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ABSTRACT 


A two color laser speckle shift strain measurement system based on the technique of Yamaguchi was 
designed. The dual wavelength light output from an Argon Ion laser was coupled into two separate 
single-mode optical fibers (patchcords). The output of the patchcords is incident on the test specimen 
(here a structural fiber). Strain on the fiber, in one direction, is produced using an Instron 4502. Shifting 
interference patterns or speckle patterns will be detected at real-time rates using 2 CCD cameras with 
image processing performed by a hardware correlator. Strain detected in fibers with diameters from 21 
microns to 143 microns is expected to be resolved to 15 fie. This system was designed to be compact 
and robust and does not require surface preparation of the structural fibers. 
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1 INTRODUCTION 


A nonintrusive strain measurement technique is required due to the small surface area of the fiber 
test specimens. In the past, extenso meters 1 or markers were used for visual observation, but these 
techniques lack the proper gage length and/or resolution to accurately determine the strain in fibers. 
A two-color simultaneous strain measurement system has been previously demonstrated 2 based on Ya- 
maguchi’s technique. 3 Laser speckles show displacement caused by deformation of a laser illuminated 
surface. Shifts in speckle patterns are directly related to the strain, translation, and rotation compo- 
nents of the deformation. A non-contact strain measurement system was designed to selectively detect 
strain in a fiber as it is pulled from one end while the other is held stationary using a driven tensile 
test machine (Instron 4502). 



2 THEORY 


The optical geometry of the system is illustrated in Figure 1. The S points mark the center of 
curvature of the incident laser wavefronts. Object deformation a(x, $/, z) of the diffuse surface causes 
speckle displacements A [Ax, Ay] at the two observation planes a distance V Q away from the object 
plane. Laser positioning lenses (not shown) provide a means of adjusting Ls , the radius of curvature. 
A far-field diffraction pattern is formed at the observation plane on the assumption that V 0 is much 
larger than the spot size of the incident beam. 



Figure 1: Geometry of optical system. 

The speckle displacements are given by 3 

Ax = -a x J^^Sx “ 1) + i* — 1 “ a y + Ixly “ &Z + Wz) 

~~L 0 [^zi(^i + l x ) + £xy{lsy 4* ly) 4- Qz{l$y 4- ly) — £ly(lsz 4* J Z )] (1) 

where ar(a x , a y , a z ) and Q(Q x ,Q y ,Q z ) are the translation and rotation vectors, respectively, and 
the components of and H are the translational and rotational components of rigid body motion, 
respectively. The unit vectors representing the direction of the center S and the center 0 are denoted 
by ls(lsxJsy, hz) and l(l x , l y , l z ), respectively. The linear x and y components of the strain tensor are 
designated e xx and e yv , respectively, and e xy is the in-plane shear component. 

The two laser sources, at angles of incidence of 9s and — 9s shown in Figure 1, are of different 
wavelengths (colors). Laser line filters F(0 5 ) and F(— 9 S ) each discriminate against the other source 
and unwanted background radiation. 
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For the angle of incidence 9s and observation angle 0 O in the x-z plane, 

l 3 = (sin 0s,O, cos 0$), and 
l = (sin0 o ,O,cos0 o ). 

The X component of A(0s,0 o ) is given by 


A(0 S} 9 o ) x = —a* (sin 0| — 1) + sin 0* — 1 j - a z ^ sin 0s cos 0s -h sin 0 O cos 0 O J 

-L' 0 [e xx (sin 0s 4* sin 0 O ) - Ct y (cos 0s + cos 0 O )] . (4) 

To eliminate the dependence of the rigid body motion component, the difference between the speckle 
displacements A(0s,0 o )x and A(— 0s,--0 o )x is calculated. That is, 


A Ax = A(0s, 9 0 )x — A{—9s, 


which after substitution becomes 


such that 


A A x = — 2L^e xx (sin 6 So 4* sin 0 o ) - 2a z \ sin 9 a cos 6 S 4- sin 0 O cos 6 0 


2 L r 0 (sin 0s 4 sin 0 o ) 


0 = j2- sin 9 S cos 0 5 4- sin 0 O cos 0 O J . (8) 

The measurement of linear strain through equation (7) is independent of the rigid body motion com- 
ponents in equation (1) as long as equation (8) is satisfied. This optical constraint can be satisfied 
experimentally or calculated through knowledge of the laser beam characteristics. From equation (7), 
the strain in the x direction is determined from the differential speckle displacement, speckle shift 
(A Ax) , and the geometry of the optical arrangement. 

Equation (7) differs from Yamaguchi’s e xx equation by the addition of the sin0 o term in the denom- 
inator. Further, in Yamaguchi’s case, equation (8) is met by having = 0, meaning the laser beam 
waist is incident on the object plane. In the two-color simultaneous system the linear arrays record the 
speckle patterns at the same time, improving the recording speed by a factor of two over a sequential 
recording system. 

It has been suggested that the linear arrays be placed one above the other rather than the present 
side by side configuration. For this proposed arrangement, the angle of incidence, 0s, is still in the x-z 
plane, but the observation angle 0 o is now in the y-z plane. Thus, l = (O,sin0 o ,cos0 o ) , and 


A (9 s ,0 o ) x = -o. x ^(sin#! — 1) — 1 - a* ^(sinfl.cosfl,) 

L-ks J 

~L' 0 [£„(sin 0 S ) + e xy (sm0 o ) + il z (sin0 o ) - Cl y (cos6 s + cos0 o )] . (9) 


Substituting equation (9) into equation (5) yields 


A.4* = -2L' 0 £ IX (sir\8s) -2a z — * -s\n6 s cos6 s - 2L' 0 [£ iy (sin 0 O ) + fi z (sin0 o )] . (10) 

L s 
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WAVELENGTH (nm) 

PERCENT OF TOTAL POWER 

457.9 

8% 

476.0 

12% 

488.0 

20% 

496.5 

18% 

501.7 

5% 

514.5 

43% 


Table 1: Percent of Argon Ion Laser Output Power at Each Wavelength 


From equation (10) it is clear that for this configuration, no value of can cancel the shear strain 
and all of the components of rigid body motion. 


A correlation algorithm is used to determine A A*. In effect, it overlays the shifted speckle pat- 
tern onto the reference pattern and varies the relative positions until the maximum product sum is 
determined. This maximum corresponds to the speckle displacement between the two patterns. This 
is represented by 


***>-;/ 


f(X)*g(X + AX)dX, 


( 11 ) 


where f(X) is the reference speckle pattern, g(X 4- AX) is the shifted speckle pattern and AX is the 
amount of the shift. The integral is computed repeatedly for increasing values of AX. The value of 
AX for which R is a maximum is the amount of shift in Ax between the before-strain and after-strain 
speckle patterns. 


This correlation is performed between the reference and shifted patterns from the 6 a incident beam, 
and then again for patterns from the — 6 S incident beam. The difference between the values of Ax is 
then used in strain equation (7) to calculate the surface strain. 


3 EXPERIMENTAL SETUP 


3.1 Optics 

A schematic of the optical system used to couple the output of an Argon Ion laser (operating in 
multi-line mode) into two single-mode fibers is shown in Figure 2. The coupling optics are mounted on 
a small optical bench which is then mounted vertically in a test cabinet. The outputs from the optical 
fibers are inserted into the sensor head having the geometry previously shown in Figure 1. 

As illustrated in Figure 2, the laser light output is reflected off a mirror onto a dispersion prism 
which separates the beam into its component wavelengths. This light is reflected off the two mirrors 
and coupled into two single- mode fibers. The power spectrum of the laser is given in Table 1, showing 
that the logical choice for specimen illumination colors are blue and green, 488 nm and 514.5 nm. An 
attenuator is placed in front of the green fiber coupler to reduce its intensity so that the blue and green 
have identical power. To eliminate undesired wavelengths coupling into each of the optical fibers, a slit 
is placed directly in front of each. 

The outputs of the optical fibers are connected to the sensor head assembly, several feet from the 
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Figure 2: Optics coupling laser light into two single-mode fibers. 


coupling optics. This assembly consists of a collimator and focusing lens and was designed for the beam 
waist to be incident on the fiber being strained. Approximately 1 W of collimated laser power must be 
coupled out of the optical fibers to illuminate the specimen. 

The two colors are incident on the test fiber at an angle of $s — 30° . Light reflected off the fiber being 
strained travels to two CCD array detectors, each with an optical filter to block unwanted wavelengths 
(See Figure 1). Because two CCD’s are required, there is an angle (0 O ) between the sample normal and 
each detector. 


3.2 Electronics 

A block diagram of the system setup is shown in Figure 3. The speckle pattern is detected by two 
linear CCD array boards controlled by the real-time hardware correlator. The correlator assembly is 
comprised of a control board and two correlator slice boards. The control board also interfaces to the 
IBM PC and stores data frames, initiates correlator sequences, and selects which stored frames will be 
used for the reference fields. Two correlator slice boards calculate the correlation for each CCD array. 
The resulting speckle displacements observed by each CCD array are converted to analog signals which 
are simultaneously digitized with the analog load and extension signals from the Instron. The control 
functions of the Instron are interfaced to the IBM PC via a GPIB card. 

The data is hard-limited 4 to take advantage of the massive parallelism of an architecture based on 
fuzzy pattern comparator ICs 5 . 6 Hard-limiting is achieved when a field of data from the CCD array 
is compared to the mean of the previous field. Data values equal to or greater than the mean are 
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Figure 3: Block diagram of electronics. 


assigned a value of one, otherwise a value of zero is assigned. While this comparison is taking place, 
the hardware is simultaneously calculating the mean for the next data set comparison. Given 8 fuzzy 
comparator ICs with 8 accumulators each, 64 bins of the correlator function are calculated as the 
Hamming distance between the incoming data set and a reference at each of the 64 displacements. The 
minimum of the Hamming function asymptotically approaches the maximum of the correlator function 
given in equation (11). The output of the correlator yields the speckle displacement, Ax- Calculating 
the difference in the displacement between the two beams provides the speckle shift, AAx - From this, 
the strain can be directly calculated using equation (7). 


4 RESULTS 

Using data acquired from a previous system setup, 7 it was shown that the minimum of the Hamming 
function asymptotically approached the maximum of the correlator function. This validated our choice 
of algorithm and hardware architecture. Furthermore, random data processed through hardware con- 
sisting of a fuzzy pattern comparator and the mean value calculating hard-limiting circuits produced 
the output shown in Figure 4. Here, the minimum corresponds to unshifted 1-bit data sets, which is 
analogous to the maximum for AX — 0 in equation (11). 
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Hamming Distance vs. Shift Displacement 



Shift Dispiac ement 


Figure 4: Hamming distance vs. shift displacement 


Comparison of Two Techniques 



Comparison of two techniques. 

Figure 5 compares actual stress vs. strain data from the original speckle rig 8 which used equation 
(11) on 8-bit speckle data to the same data hard-limited and processed by the Hamming distance algo- 


7 




rithm. Results shown agreed, indicating that the algorithm cho6en can be suitably substituted for the 
correlator -based system. Using the traditional correlator algorithm requires a great deal of computa- 
tional overhead which results in a time and/or hardware penalty. The 1-bit Hamming algorithm uses 
smaller data paths, less computational intensity, and allows for the use of commercial off the shelf ICs. 
Thus, this method efficiently processes the data to exceed the real-time control rate of the Instron, 
thereby allowing for off-line characterization of samples. 


5 SUMMARY 


A two-color speckle strain system has been designed. The optics have been completed and tested. 
Proof of concept electronics have been tested and a system designed to achieve real-time data rates. In 
the future, an oven will be placed around the fiber mount of the Instron such that strain measurements 
can be taken at temperatures exceeding 1000°C degrees. 
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